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I.  Introduction 

A  computationally  tractable  and  physically  clear  model  for 
rotational  autoionization  of  dipole-bound  states  of  molecular  anions 
has  recently  been  put  forth^^  by  Clary.  In  this  model,  the  "extra" 
electron,  originally  bound  and  locked  to  a  molecule-fixed  axis  system, 
gains  energy  from  the  molecule’s  rotational  degrees  of  freedom  and 
is  eventually  ejected,  leaving  the  molecule  with  less  rotational 
energy.  The  energy  levels  of  the  anion  relative  to  those  of  the  neutral 
molecule  and  the  autoionization  lifetimes  of  these  levels  are 
predicted  within  Clary’s  model  in  terms  of  the  dipole  moment  and 
rotational  constants  of  the  neutral  molecule.  Although  this  model  has 
been  remarkably  successful  in  interpreting  trends  in  the 
experimental  spectra^^  of  these  metastable  states,  certain  aspects  of 
the  data  require  enhancement  of  the  original  model.  It  is  in  this 
spirit  that  the  developments  presented  here  are  put  forth. 

II.  Overview  of  Clary's  Original  Rotationally  Adiabatic  Model 

A.  The  Angular  Basis  Functions  and  the  Rotationally  Adiabatic  Curves 
To  introduce  the  main  concepts  and  notation  of  the  rotationally 
adiabatic  (RA)  modeU^,  it  is  adequate  to  consider  an  electron  bound 
to  a  rotating  diatomic  species^  whose  vibrationally  averaged 
rotational  constant  is  denoted  by.  To  describe  states  of  a  given  total 
angular  momentum  J  and  projection  M  along  the  lab-fixed  z-axis,  the 
following  basis  of  angular  functions  is  formed 

I  J,M;  j,l  >  =  Zm  <  j,m,l,M-m  I  J,M  >  I  j,m  >  I  l,M-m  >. 

These  functions  consist  of  coupled  products  of  molecular  rotor 
(I  j,m  >)  and  electronic  (I  l,M-m  >)  angular  functions.  A  model 
Hamiltonian  containing  a  molecular  rotor  part  (by  j^),  an  electron- 
molecule  potential  V,  and  an  electronic  angular-plus-radial  kinetic 
energy  part  is  also  introduced 

H  =  by  j2  +  V  +  12  /2meR2  -h2/2me  1/R2  (d/BR  (  R2  d/dR)). 

The  electron-molecule  potential  is  tal  n,  for  these  dipole-bound 
anion  states,  to  be  of  the  form 
V  (R,0)  =  p  cos(0)/R2  F(R), 

where  F(R)  includes  the  effect  of  the  electron-polarizability 
interaction  and  a  "cut  off"  function  designed  to  smoothly  remove  the 
singularity  in  the  potential  as  R  approaches  zero. 

The  key  construct  of  the  original  RA  model  as  put  forth  by 
Clary  is  the  set  of  so-called  rotationally  adiabatic  potential  curves 
generated  by  diagonalizing,  within  the  angular  basis  described  above, 
the  Hamiltonian  in  the  absence  of  electronic  radial  motion 
HO  =  by  j2  V  -h  12  /2me 
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for  fixed  electron-molecule  distance  R.  For  each  fixed  J  value,  Clary 
found  that  the  lowest  rotationally  adiabatic  potential  curve 
generated  in  this  manner  was  attractive  at  small  R,  had  a  "barrier"  at 
intermediate  R  values,  and  approached  the  energy  of  the  j  =  0  rotor 
at  large  R.  The  physical  origins  of  these  attributes  of  this  lowest 
potential  are  discussed  below  in  Sec.  III. 

B.  The  Dipole-Bound  State  Energies  and  Lifetimes 

The  RA  model  evaluates  the  dipole-bound  anion's  electronic 
energies  and  autoionization  lifetimes  by  using  the  lowest  RA 
potential  curve  (call  it  Eq  (R))  to  solve  for  the  radial  motion  of  the 
electron.  The  energies  of  the  resultant  quasi-stationary  states  and 
their  lifetimes  (which  are  governed  by  the  rate  of  tunneling  through 
the  barrier  in  Eq  (R))  are  obtained  from  this  radial  Schrddinger 
equation's  solution.  Much  of  the  solution  was  carried  through 
analytically  la  by  Clary,  and  many  of  the  trends  in  decay  rates  as 
functions  of  rotational  quantum  numbers,  dipole  moment,  and 
isotopic  substitution  were  beautifully  explained  by  this  model. 

III.  Observations  Leading  to  Modification  of  the  Clary  Model 
A.  The  Nature  of  Rotationally  Adiabatic  Energies 

It  is  important  to  understand  the  physical  nature  of  the  lowest 
rotationally  adiabatic  potential  Eo(R)  and  of  the  corresponding 
angular  eigenfunction  <l>o(R)  formed  as  a  linear  combination  of  the 
above  1J,M;  j,l  >  basis  states  with  coefficients  that  depend  on  R: 

<I)o(R)  =  2:j,iCo.j,i  (R)  lJ,M;j,l>. 

This  angular  function,  multiplied  by  the  solution  (F(R))  of  the  radial 
Schrodinger  equation,  forms  the  full  electronic-rotational 
wavefunction  within  this  model:  'F  =  Oq  F(R).  To  understand  the 
nature  of  the  couplings  among  the  IJ,M;  j,l  >  that  give  rise  to  Eo(R)  and 
<l>o(R),  consider  the  "diabatic  energies"  given  by  the  expectation 
values  of  HO  for  each  of  the  angular  basis  states  IJ,M;  j,l  >  : 

E0j.i  (R)  =  <  J,M;  j,ll  HO  IJ,M;  j,l  >  =  by  jU+D  +  1(1+1) ■h2/2meR2  . 
Because  the  angular  dependence  of  V  is  of  the  form  cos(0),  the 
average  value  of  V  taken  over  any  single  IJ,M;  j,l  >  function  vanishes. 
Thus,  in  the  absence  of  couplings  among  the  angular  basis  states,  one 
has  a  picture  of  diabatic  potentials  EOj  i  (R)  which  undergo  multiple 
crossings  as  depicted  below  in  Fig.  1. 

Couplings  among  angular  basis  states  which  have  the  same  J 
and  M  values  (  these  two  quantum  numbers  are  conserved  by  the 
full  Hamiltonian  H)  but  which  differ  by  unity  in  their  j  (and  hence  1) 
quantum  numbers  can  be  realized  through  the  angular  dependence 
of  V(R,0).  At  small  R,  where  EOj^q  and  EOj.ij  undergo  a  crossing,  the 
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two  basis  states  with  j=J,  1=0  and  j=J-l,l=l  mix  to  form  a  hybrid  state 
Oo  (i.e.,  they  combine  to  dominate  the  sum  Z  j,i  Co.j.i  (R)  I  J,M;  j,l  >  at 
such  values  of  R)  whose  average  interaction  with  the  electron 
<<1)0  IVI  Oo>  is  non-zero.  This  hybrid  state  consists  of  an  electron  in  an 
1=0,  1=1  hybrid  orbital^  bound  and  locked  to  the  diatomic  molecule’s 
axis  with  the  molecule  rotating  in  a  manner  to  produce  total  angular 
momentum  J.  Binding  is  caused  largely  by  attraction  of  the  electron 
to  the  positive  end  of  the  molecule’s  dipole  field. 

At  larger  values  of  R,  and  E^j.2,2  undergo  a  crossing; 

here  ^0  picks  up  a  strong  component  of  j=J-2,l=2  character,  and  Eo(R) 
evolves  to  track  EOj-2,2  (R)  beyond  this  crossing  until  the  next 
crossing  is  reached.  Subsequent  crossings  between  EOj.n.n  and 
E®j-n-l,n+l  occur  at  Rj-n,J-n-l  =  {(n+l)/(2b(J-n))  at  energies  of 

Ej-n.J-n-i=b(J-n)(J+l)).  These  crossing  energies  lie  above  the  neutral’s 
rotational  state  J-n  (whose  energy  is  b(J-n)(J-n-i-l))  by  an  amount‘d 
nb(J-n)).  As  subsequent  crossings  are  reached,  couplings  between  the 
j=J-n,  l=n  component  of  <I)o,  which  exists  as  dominant  from  the 
preceeding  crossing,  and  the  j=J-n-l,  l=n+l  component  give  rise  to  a 
new  (j=J-n-l,  l=n+l)  contribution  to  <I>o  and  to  a  further  lowering  of 
Eo  (R)  as  this  potential  moves  from  tracking  E^j.^.n  to  track 
E^J-n-l.n+l  • 

B.  Problem  with  Using  Eq  (R)  at  Large  R 

One  difficulty  of  the  original  rotationally  adiabatic  model 
should  now  be  clear;  by  forming  Eo(R)  in  this  manner  and  solving  for 
the  radial  motion  of  the  electron  on  the  single  potential  surface 
Eo(R),  one  restricts  the  autoionization  process  to  yield  a  state 
^=00  F(R)  which  contains  only  j=0,  1=J  components  at  large  R  because, 
<1)0  contains,  at  any  R,  primarily  that  j,l  component  arising  from  the 
last  crossing  encountered.  This  model  does  not  allow  for  the 
formation  of  the  numerous  final  states  (e.g.,  j=J-l,  3-2,  ...3,  2,  1,  0  may 
be  energetically  accessible)  that  may  be  "open"  (i.e.,  accessible).  To 
extend  the  model  in  a  manner  that  will  allow  final-state  branching 
ratios  to  be  addressed  requires  treatment  of  the  multi-channel 
nature  of  the  electron-ejection  dynamics. 

It  is  natural  to  wonder  how  this  model  was  able  to  provide  so 
much  insight  into  the  energy  level  patterns  and  (total)  decay 
lifitimes  of  the  dipole-bound  states  which  it  is  designed  to  mimic. 

The  answer  lies  in  the  observation  that  the  positions  (i.e.,  energies) 
and  lifetimes  of  the  states  obtained  by  solving  the  radial  Schrodinger 
equtation  on  Eo(R)  depend  primarily  on  the  form  of  Eo(R)  where  this 
potential  is  attractive  (for  the  positions)  and  on  the  height  and  width 
of  the  barrier  at  these  energies(for  the  lifetimes).  Beyond  the  barrier. 
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where  the  electron  no  longer  is  bound,  is  it  unreasonable  to  make  the 
rotationally  adiabatic  approximation.  Thus,  the  form  of  and  use  of 
Eo(R)  is  correct  at  those  R  values  (i.e.,  small  R  and  throughout  most  of 
the  barrier  region)  that  affect  the  positions  and  lifetimes  of  the 
states.  The  fact  that  Clary's  rotationally  adiabatic  Eo(R)  does  not 
provide  an  adequate  description  of  the  dynamics  beyond  the  barrier, 
where  a  multi-open-channel  picture  is  needed  and  where  the 
electron  is  no  longer  locked  to  the  molecular  axis,  has  very  little 
influence  on  the  model's  predicted  energies  and  total(  not  branching- 
specific)  lifetimes. 

C.  The  Predicted  Rotational  Energy  Level  Patterns 

In  addition  to  the  need  to  augment  the  original  Clary  model  by 
including  a  multi-open-channel  treatment  at  large  R,  the  rotational 
energy  level  pattern  obtained  from  the  original  model  is  not  accurate 
enough  to  characterize  certain  aspects  of  the  experimentally 
observed  lifetime  trends.  In  particular,  the  original  model  expresses 
the  rotational  energies  of  the  dipole-bound  states  in  terms  of  only 
one  set  of  rotational  constants-those  of  the  daughter  neutral 
molecule.  The  anion's  energies  are  expressed  in  terms  of  these 
rotational  constants  and  the  dipole  moment  of  the  neutral  molecule. 
As  given  explicitly  in  Eq.(14)  of  ref. la,  the  energy  levels  of  the 
dipole-bound  anion  states  (e.g.,  Eres(K’,J')  =  -C(b)  +  bJ'(J’+l)  +  (a-b)  K'2, 
for  a  symmetric  top  case)  are  related  to  those  of  the  neutral  by  a 
single  negative  constant  -C(b)  which  is  independent  of  rotational 
quantum  numbers.  This  constant  shift  (i.e.,  the  electron  affinity  EA) 
in  the  anion  energies  relative  to  those  of  the  neutral  is  an  integral 
part  of  the  original  model.  This  gives  rise  to  a  situation  in  which: 

1.  For  an  electron  affinity  EA  of  the  dipole-bound  state,  the  lowest 
state  of  the  anion  which  can  detach  has  J'  (J’+l)  b  =  EA.  This  state  can 
detach  only  via  a  AJ  =  J'  process. 

2.  As  one  moves  to  higher  J  values  within  the  dipole-bound  anion's 
rotational  "ladder",  smaller  and  smaller  AJ  processes  become  open 
(see  the  table  below). 

3.  Eventually,  for  high  enough  J,  one  reaches  a  situation  in  which  AJ  = 
1  is  open,  but  one  never  realizes  the  AJ  =  0  situation.  From  these 

J  values  on,  AJ  =  1  remains  open. 

Such  a  progression  from  initially  large  AJ  through  smaller  AJ 
values  as  J’  increases  is  often  seen  experimentally;  however,  there 
are  interesting  cases^  in  which  the  pattern  of  AJ  vaules  is 
qualitatively  different.  It  turns  out  that  substantial  errors  in 
calculated  detachment  rates  can  occur  when  the  original  model 
incorrectly  predicts  the  anion-neutral  state  spacings  in  regions  where 
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AJ  is  changing.  For  example,  for  J'  values  where  AJ  =2  transitions  are 
experimentally  seen  to  be  closed  but  nearly  open,  the  original  model 
may  (incorrectly)  predict  AJ  =  2  to  be  open.  It  is  important  to  correct 
such  small  errors  in  the  state  energies  predicted  by  the  original 
model  in  these  special  cases  where  AJ  is  incorrect  because  these 
errors  may  lead  to  qualitatively  incorrect  neutral-state  branching 
ratios.  As  shown  below,  these  same  corrections,  combined  with 
modification  of  the  model  introduced  to  more  properly  treat  the 
asymptotic  multi-channel  issues,  allow  one  to  better  understand  the 
sudden  rises  in  detachment  rate  that  can  accompany  new  channel 
openings  (i.e.,  smaller  AJ  values  becoming  open). 

IV.  The  Modified  Rotationally  Adiabatic  Model 

To  overcome  the  lack  of  multiple-open-channels  at  large  R  and 
to  represent  the  anion-neutral  energy  level  patterns  in  a  manner 
which  allows  progressive  channel  openings  to  be  consistent  with 
different  values  for  b^  and  b*,  the  original  rotationally  adiabatic 
model  can  be  modified  as  follows: 

A.  Beyond  the  crossing  point  Rj.n,j.n-l  =  {(n+l)/(2b(J-n))  of  the 
first  open  channel  (J-n),  the  full  multitude  of  open-channel  surfaces 
are  brought  into  consideration  and  the  electronic  function  is 
permitted  to  develop  amplitudes  for  each  such  channel.  In  the 
simplest  treatment  of  this  aspect  of  the  problem,  a  Landau-Zener 
model  can  be  used  to  estimate  the  amplitudes  of  the  final-states 
whose  crossing  energies  Ej-n,j-n-l=b(J-n)(J-i-l))  occur  below  the 
dipole-bound  state's  energy  E.  Decay  into  final  states  whose  crossing 
energies  lie  above  E  can  be  treated  perturbatively  as  described 
below. 

B.  To  more  accurately  represent  the  relative  orderings  of  the  dipole- 
bound-anion  and  neutral-molecule  rotational  states  and  thereby 
achieve  a  more  correct  description  of  the  AJ  pattern,  the  energies 
obtained  by  solving  the  radial  Schrodinger  equation  on  Eo(R)  are 
"fine  tuned"  by  simply  correcting  them  to  fit  the  experimentally 
observed  pattern.  The  necessary  corrections  are  small  and  should  not 
substantially  affect  the  tunnelling  rates  calculated  within  this  model 
because  the  width  of  the  barrier  is  only  weakly  dependent  on  the 
position  (energy)  of  the  state  for  such  small  energy  variations.  The 
primary  effect  of  such  energy  shifting  is  to  close  or  open  channels 
which  the  original  model  (with  b®  =  b')  incorrectly  identified.  These 
channel-opening  corrections  can  have  substantial  affects  on  the 
calculated  detachment  rates  and  branching  ratios. 
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V.  Implementation  of  the  Modified  Model 

After  "fine  tuning"  thi  rotational  levels  of  the  dipole-bound 
anion  predicted  by  the  original  model  (labeled  by  J  or  J,K  for 
diatomic  or  symmetric  top  species,  respectively)  so  as  to  correctly 
identify  the  first  energetically  available  neutral-molecule  state  (label 
it  Ji  =  J-n)  for  each  anion  state  J,  the  total  rate  of  autoionization  is 
first  computed  as  in  the  original  rotationally  adiabatic  model. 
Beginning  with  the  state  J-n,  branching  ratios  are  then  calculated.  For 
open  channels  whose  crossing  energies  lie  below  the  anion's  energy 
E,  the  Landau-Zener  (LZ)  method  as  detailed  below  is  used.  For  open 
channels  whose  crossing  energies  Ej-n.j-n-l  lie  above  E  (such 
channels  are  denoted  "frustrated"),  the  perturbative  method  given 
below  is  used.  This  treatment  of  the  frustrated  channels  also 
provides  a  correction  to  the  overall  tunnelling  detachment  rate 
predicted  by  the  original  model,  arising  from  channel-opening 
effects. 

A.  Perturbative  Treatment  of  Frustrated  Channels 

For  open  channels  whose  crossing  energies  Ej.n,j.n-l  lie  above 
E,  one  can  estimate  the  branching  ratios  between  the  corresponding 
pairs  of  states  (j=J-n,  l=n  and  j=J-n'n,  l=n-»-)  by  considering  the  2  x  2 
couplings  between  successive  pairs  of  states.  As  the  diabatic  surfaces 
and  .n+l(R)  couple  through  Vj.n,j-n-l(R).  the  mixing 

coefficients  Cj-n  and  Cj-n-i  corresponding  to  the  two  coupled  states 
will  vary  with  R.  For  the  highest  energy  frustrated  state,  the 
wavefunction  Oo(R)  brings  to  the  Rj-n.J-n-l  crossing  a  dominant 
amplitude  for  the  j=J-n  state.  To  estimate  the  amplitude  of  the 
j=J-n-l  state  once  this  crossing  point  has  been  passed,  perturbation 
theory  can  be  used  and  gives: 

Cj-n-l/Cj-n  =  Vj-n.J-n-l(R)/[E-EOj-n-l(R)  ]. 

The  probability  that  the  system  would  "hop"  from  the  diabatic  state 
with  j=J-n  to  that  with  j=J-n-l  as  it  passes  through  this  crossing  is: 

Pn.n+l  =  I  Cj-n-1  \'^l  [I  Cj-n-l  |2  Cj.„  |2]; 
the  probability  that  it  would  remain  on  the  state  with  j=J-n  is 
(1  ■  Pn,n+)- 

The  energy  can  reasonably  be  expressed  as  E  =  b(J-n)(J-n-i-l) 
+fnb(J-n)  since,  by  assumption,  this  dipole  bound  state  lies  above  the 
j=J-n  state  of  the  neutral  and  below  the  crossing  which  occurs  at 
Rj-n.j-n-l  =  {  (n+l)/(2b(J-n))  }  1/2  with  energy  Ej-nj-n-l  =  b(J-n)(J-t-l). 
Here  f  parameterizes  the  state  energy  E  in  terms  of  the  spacing 
(nb(J-n))  between  Ej-n,j-n-l  and  Ej.n.  Just  beyond  the  crossing  point, 
where  the  wavefunction’s  flux  bifurcates,  the  square  of  the 
amplitude  ratio  can  thus  be  estimated  as: 
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I  Cj-n-i  |2/  I  Cj-n  |2  =  4(V0j.„,j.„.i)2/[(i .f)2(n+ 1  )2n2] , 
where  the  matrix  element  of  V  between  the  two  diabatic  states 
Vj-n,J-n-l  has  been  written  in  as  V0/R2j.n,j.n-i.  These  estimates  for 
the  amplitude  ratios  lead  to  the  following  predictions  for  the  hopping 
probabilities  involving  frustrated  states: 

Pn,n+1  =  4(V0j.„j.„.i)2/[n2  (n+l)2(l-f)2  +  4(V0j.n,j.„.i)2],  and 

1-Pn,n+1  =n2  (n+l)2  (l-f)2/  [n2  (n+l)2(l-f)2  +  4(V0j.n,j.n-l)2] . 
After  computing  the  probabilities  for  the  highest-energy  frustrated 
channel,  one  moves  to  the  next  frustated  channel  (if  there  is  one)  and 
computes  Pn+i,n+2  in  the  same  manner. 

As  stated  above,  flux  that  does  not  "hop"  to  one  lower  j -value 
at  a  given  crossing  remains  on  the  higher  diabatic  curve.  Once  on  this 
curve,  its  amplitude  must  still  tunnel  through  the  remaining 
classically  forbidden  region;  this  tunnelling  further  reduces  the 
probability  of  eventually  exiting  in  this  product  channel  by  an 
amount  exp(-2KL).  Here,  L  is  the  "extra  distance"  that  the  electron 
must  tunnel  from  Rj-n,J-n-l  to  where  it  exits  the  classically  forbidden 
region  on  the  E0j.n,n(R)  diabatic  potential.  The  latter  distance  can  be 
calculated  by  first  expressing  the  dipole-bound  state  energy  E  as 
above:E  =  b(J-n)(J-n+l)  +  fnb(J-n).  Here,  b(J-n)(J-n-hl)  is  the 
asymptotic  energy  of  the  j=J-n  frustrated-  channel  state  and  nb(J-n) 
is  the  energy  gap  between  this  asymptote  and  the  Ej-n.j-n-l  crossing 
energy.  Thus,  f  characterizes  the  fraction  of  the  way  between  the 
asymptote  and  the  frustrated-crossing  energy.  Setting  this  energy 
equal  to  b(J-n)(J-n+l)  -i-  n(n+l)/2R2  yields  the  value  of  R  at  which 
flux  exits  the  classically  forbidden  energy  on  E^j.n.nCR): 

Rexit  =  {(n+l)/(2b(J-n)f}l/2  . 

The  K  value  is  given  by 

K  =  {2m[  Ebarricr  *  E]}^/2  ^ 

where  Ebarrier  is  the  potential  through  which  the  electron  is 
tunnelling  between  Rj-n.j-n-l  and  Rexit  •  The  resultant  tunnelling 
correction  exp(-2icL)  is  then  calculated  as: 

Xn(f)  =  exp{-2(n(n-»-l))l/2  [((l-0/f)l/2  -  (l-f)l/2]}. 

Clearly,  as  the  fraction  f  approaches  zero  (i.e.,  when  the  channel  is 
barely  open),  Xn  approaches  zero;  little  flux  is  able  to  tunnel  through 
the  E0j-n,n(R)  potential  to  yield  products  in  the  j=J-n  channel.  In 
contrast,  when  f  approaches  unity,  Xn  approaches  unity  and  all  the 
flux  that  remains  on  the  j=J-n  diabatic  surface  yields  product  in  this 
channel. 

As  shown  below,  the  product  yield  for  the  channel  with  j=J-n  is 
(l-Pn,n+l)Xn  .  For  the  frustrated  channel  with  j=J-n-l,  the  yield  is 
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Pj-n.J-n-l(l-  Pj-n-i.J-n-2)Xn+i  '  the  fraction  that  hopped  to  j-J-n-1  at 
the  first  crossing  and  subsequently  remained  on  j=J-n-l  at  the 
second  frustrated  crossing.  The  fraction  Pj-n,J-n-l  Pj-n-l,J-n-2  is 
predicted  to  have  exited  the  second  crossing  on  j=J-n-2.  This 
recursive  process  is  continued  until  the  first  open,  non-frustrated 
channel  is  reached  (i.e.,  the  first  channel  j=J-n*  for  which  E  exceeds 

B.  Open  Non -Frustrated  Channels 

Once  a  channel  for  which  E  lies  above  the  crossing  energy  is 
encountered,  the  Landau-Zener  (LZ)  method  can  be  employed  to 
calculate  subsequent  hopping  probabilities.  The  fact  that  the 
potential  V(R,0)  depends  on  0  as  cos(0)  implies  that  flux  with  j=J-n  , 
l=n  can  "hop"  to  the  j=J-n-l,  l=n  +1  surface  or  remain  on  the  j=J-n  , 
l=n  surface.  That  fraction  of  the  flux  which  remains  on  the  j=J-n  ,  l=n 
surface  must  remain  on  this  surface  all  the  way  to  R  =  «»  because  the 
potential  V(R,0)  can  not  couple  j=J-n  ,  l=n  to  j=J-n-2,  l=n+2  or  to  any 
other  diabatic  surface  which  it  encounters  via  crossings  at  larger  R 
values.  That  fraction  of  the  flux  that  "hopped"  to  the  j=J-n-l,  l=n+l 
surface  can,  at  the  next  crossing  of  diabatic  surfaces,  either  remain 
on  this  surface  (again,  all  the  way  to  R=«»  because  the  potential  can 
no  longer  couple  j=J-n-l,  l=n  +1  to  other  diabatic  states)  or  hop  to  the 
j=J-n-2,  l=n  +2  surface.  This  progression  of  surface  hoppings  proceeds 
all  the  way  to  the  j=l,  1=J  -1  and  j=0,  1=  J  crossing.  At  successive 
crossings,  the  couplings  grow  weaker  because  of  the  l/R^ 
dependence  of  V(R,0). 

Within  the  Landau-Zener  method^,  the  probability  (Pn,n+l)  for 
"hopping"  from  the  diabatic  surface  with  j=J-n  ,  1=  n  to  that  with  j=J  - 
n-1,  1=  n+1  is  given  in  terms  of  matrix  elements  (Vj-n.j-n-l)  of  the 
potential  V(R,0)  between  the  two  appropriate  angular  basis  states, 
the  difference  in  "slopes"  (ASn,n+l)  of  the  two  crossing  diabatic 
surfaces  at  the  crossing  point  (Rj-n,J-n-l)  and  the  classical  radial 
velocity  of  the  electron  (  dR/dtn.n+l)  at  the  crossing: 

Pn.n+1  =  l-exp(-Q),  where 
Q  =  I  Vj-n.j-n-l  |2/{h  (dR/dtn.n+l)  I  ASn.n+1  1). 

Within  the  original  Clary  model,  the  crossing  points  and  slope 
differences  can  be  computed  and  are  given  by: 

Rj-aJ-n.l  =  ((n+l)/2b(J-n))l«  and 
AS„.„+i  =.  2(n+l)  (2b(J-n)/(n+l))3«  . 

The  potential  coupling  matrix  elements  Vj-nj-n-l  are  given  by  Clary 
as*®; 


Vn.n+1  =(n/R^){(2(J-n)+l  (2(J-n-l)+l)(2n+l)(2(n+l)+l)}'^^ 

(-l)’''*Yn  n+1  l\/-n  J-n-1  1  \  Jj-n-l  n+1  J 

\  0  0  (7  \0  0  0  /  n  J-n  1 

k 

where  the  3-j  and  6-j  symbols  arise  from  coupling  the  angular 
momenta  of  the  initial  and  final  states  through  that  carried  by  the 
cos(0)  component  of  V(R,0).  These  angular  momentum  factors  can  be 
expressed  explicitly  and  give: 

Vn.n+l  =0.20pDebyes  {(n+l)(J-n)/((n+3/2)(J-n+l/2))}  1/2/R2. 

The  radial  velocity  at  the  crossing  point  Rj-n,j-n-l  is  given  in  terms  of 
the  total  energy  E  of  the  detaching  state  and  the  functional  form  of 
the  diabatic  curve  E^j-n  n(R)  ^ 

(dR/dtn.n+l)  =  {2['  E  -  b(J-n)(J-n+l)  -  n(n+l)/2R2  ])U2 
=  {2[  E  -  b(J-n)(J+l)  ]}I/2. 

As  noted  above,  only  those  final  states  j=J-n  with  b(J-n)(J+l)  >  E  can 
be  examined  using  the  LZ  approach  as  outlined  here.  The  hopping 
probabilities  involving  frustrated  open  channels  are  computed  using 
the  perturbative  approach  detailed  earlier. 

C.  Product-Channel  Yields 

The  probability  that  the  system  hops  to  the  diabatic 
curve  with  one  lower  j  value  is  Pn,n+l  ;  the  probability  that  it 
remains  on  the  incoming  diabatic  curve  (and  thus  proceeds  all  the 
way  to  R=o«’  on  this  curve)  is  1-Pn.n+l-  Thus,  given  knowledge  of  the 
highest-energy  open  channel  (call  it  Ji),  one  evaluates  Pn,n+1  for 
n=J-Ji  .  If  Ji  is  frustrated,  then  the  perturbative  equations  are  used; 
else,  the  LZ  expressions  are  employed.  The  total  yield  of  neutrals  in 
rotational  state  Ji  is  given  by  thisl-Pn.n+l  (multiplied  by  Xj.ji  if  this 
channel  is  frustrated).  For  the  Pn,n+1  fraction  that  hop  to  the  j=  J-n-1 
diabatic  curve  (i.e.,  to  the  diabatic  curve  with  one  lower  j  value),  one 
computes  Pn+l,n+2  .  and  evaluates  the  yield  of  neutrals  in  Ji-1  as 
Pn,n+1  (l-Pn+l,n+2)  (  multiplied  by  Xj-ji-i  if  this  channel  is  also 
frustrated).  For  the  Pn,n+l  Pn+l,n+2  fraction  that  hopped  to  the  n-f-2 
curve,  one  computes  Pn+2,n+3  and  evaluates  the  yield  of  neutrals  in 
Ji-2  as  Pn,n+1  Pn+l.n+2  (l-Pn+2,n+3)  (multiplied  by  Xj.ji-2  if  this 
channel  is  also  frustrated)  .  This  recursive  process  is  continued  until 
one  reaches  the  last  crossing.  The  general  formula  for  production  of 
neutrals  in  level  Jj-k  =  J-n-  k  is: 

Y(Ji-k)  =  {nis.i,k  Pn+l-l,n+l}  (1*  Pn+k,n+k+l)Xj.ji-i(;., 
where  the  Xj.ji-k  are  the  tunnelling  corrections  that  apply  only  to 
frustrated  open  channels  that  occur.  This  procedure  gives  the 
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fraction  of  detaching  anions  that  begin  in  dipole-bound  state  J  that 
end  up  in  neutral-molecule  state  Ji-k,  where  k  ranges  from  0  through 
Ji  (when  k=Ji,  the  yield  formula  is  {ni=i,k  Pn+l-l.n+l)).  with  Ji  being 
the  first  (frustrated  or  not)  open  channel. 

VI.  Rate  Enhancements  Near  Channel  Openings 

As  one  passes  through  spectral  regions  where  new  AJ  values 
open  (e.g.,  for  several  J  values,  AJ=2  may  be  open  but  AJ  =1  is  closed; 
then,  at  higher  J,  AJ=1  opens),  the  rate  of  detachment  is  often 
observed  to  undergo  significant  increases.  The  original  rotationally 
adiabatic  model  is  able  to  reproduce  many  of  the  trends  in  lifetimes 
in  terms  of  the  dependence  of  Eo(R)  (in  particular,  the  tunnelling 
barrier  width  and  height)  on  rotational  quantum  numbers.  However, 
there  seem  to  be  effects  that  arise  in  the  neighborhoods  of  channel 
openings  which  require  extension  of  the  model. 

To  understand  these  effects,  consider  a  situation  in  which  the 
dipole-bound  state  energy  Ej.i  lies  just  below  the  J-n-1  state  of  the 
neutral  but  the  anion's  state  Ej  lies  just  above  the  J-n  state  of  the 
neutral;  in  this  case,  AJ  changes  from  n+1  to  n  as  the  anion  moves 
from  Ej.i  to  Ej  .  The  qualitatively  important  effect  which  the  original 
rotationally  adiabatic  model  neglects  has  to  do  with  the  fact  that  the 
neutral  state  J-n  is  a  "frustrated"  state  when  the  dipole-bound  state’s 
quantum  number  is  J  (in  fact,  J-n-1  and  others  may  also  be 
frustrated  if  their  crossing  energies  Ej.n-i,j-n-2  .  etc.  lie  higher  than 
the  dipole-bound  state's  energy).  In  moving  from  a  closed  channel  to 
a  frustrated  channel,  the  original  model  makes  no  qualitative  change 
in  its  ionization-rate  predictions,  which  remain  goverened  purely  by 
the  rate  of  tunnelling  through  Eo(R).  The  original  model  has  the  open 
but  frustrated  channels  contribute  no  additional  effect  to  the  total 
detachment  rate  (and  gain  zero  branching  yield,  as  a  result)  because 
no  tunnelling  occurs  to  R  values  beyond  their  crossing  points 
Rj-n.J-n-i,  etc. 

Within  the  frustrated-channel  picture,  amplitude  on  the 
neutral's  frustrated  states  are  present  in  amounts  which  are 
estimated  by  perturbation  methods.  This  analysis  also  gives  an 
estimate  of  the  fractional  increase  in  the  detachment  rate  to  be 
expected  a  one  experiences  a  new  channel  opening.  In  the  above 
example  j-.  prior  to  the  channel  opening  when  the  dipole-bound 
state's  qiw..um  number  is  J-1,  the  state  lying  just  above  E  (i.e.,  the 
j=J-n-l  state  which  AJ  =  n  in  the  above  example)  is  closed;  only 
the  frac  ion  Pn+l,n+2  in  the  j=J-n-2  state  continues  on  to  eventually 
detach  in  J-n-2,  J-n-3,  etc.  After  the  opening  when  the  dipole-bound 


state  quantum  number  is  J,  both  the  corresponding  fraction  (Pn,n+l) 
and  the  fraction  (1-  Pn,n+i)Xn  remaining  on  the  new  open  channel 
j=J-n  eventually  detach.  Thus,  the  rate  ratio  is  given  by; 

ratej/ratej.i  =  [(1-  Pn,n+l)Xn  +  (Pn,n+l)]/Pn+l,n+2- 

The  original  rotationally  adiabatic  model  ignores  the  fraction 
(1-  Pn.n+l)Xn  that  detach  via  the  frustrated  channel  and  counts  only 
the  fraction  Pn,n+l  that  detach  after  hopping  to  j=J-n-l.  For  the 
dipole-bound  state  Ej  in  which  j=J-n  is  frustrated,  the  predictions  of 
the  original  model  must  therefore  be  scaled  by  the  ratio: 

ratej/ratej, adiabatic  =  [(!■  Pn,n+l)Xn  +  (Pn,n+l)]/Pn,n+l 

=  l  +  {n2  (n+l)2  (l-02/4(V0j.„.j.n.O2} 

exp{-2(n(n+l))l/2  [  ((l-f)/f)l/2  -  (1-01/2  ]  } 

=  l  +  {n2(n+l)(l-f)2(n+3/2)(J-n+l/2)/((J-n)p2i3^bycs)) 

exp{-2(n(n+l))l/2  [((l-O/Ol/^  -  (l-O^/^J). 
where  the  expressions  given  earlier  for  the  hopping  probabilities 
among  frustrated  channels  and  for  the  V0j-n,j-n-l  have  been 
inserted.  In  fact,  if  lower  open-channel  states  are  also  frustrated 
(this  can  occur  if  their  crossing  energies  lie  above  E),  the  rate  sh>  d 
be  scaled  by  a  product  of  factors  as  given  above  where  n  is  replaced 
by  n+k  and  k  is  allowed  to  run  over  all  states  j=J-n-k  that  are 
frustrated  at  energy  E.  These  scalings  are  necessary  because  the 
original  adiabatic  model  allows  only  Pn+k.n+k+l  of  the  amplitude 
entering  each  of  the  frustrated  crossing  points  to  detach;  in  reality, 
all  of  the  amplitude  may  detach. 

VII.  Example  Application 

To  illustrate  the  effects  of  the  frustrated-channel  corrections 
described  here,  consider  dipole-bound  states  arising  from  a  molecule 
with  p  =  3.6  Debyes  (this  is  approximately  correct  for^b  CH2CN).  The 
exponential  tunnelling  component  of  the  correction  factor  varies  with 
n  and  f  as  shown  in  Table  2  below.  The  factor  W(n,J)/p2£,gjjyg5 
=6.465{n2(n4-l)(n-t-3/2)(J-n-(-l/2)/((J-n)p2Qgjjyg5))  is  given,  for  various 
n,  as  in  Table  3.  The  overall  ratio  of  the  predicted  to  rotationally 
adiabatic  rates  is  then  given  as  l+T(n,f)W(n,J)/p2£,gjjygg. 

It  is  useful  to  tabulate  quantities  in  the  above  manner  because, 
for  any  given  dipole  moment  poebycs.  one  can  easily  compute  the 
channel-opening  corrections  to  the  RA  model  in  terms  of  the 
experimentally  observed  anion-neutral  state  energy  gap  (which 
determines  f).  For  example,  consider  the  AK=0  detachment  rates 
associated  with  the  J=  31-40,  K=0  levels  of  the  dipole-bound  state  of 
CH2CN-.  These  states  all  decaylb  to  neutral-molecule  rotational  levels 
with  AJ^3.  The  anion  and  neutral  (with  AJ=3)  energy  levels  for  this 
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system,  along  with  the  corresponding  f-values,  are  given  below  in 
Table  4  (here  n=3). 

Using  a  dipole  moment^h  of  pDebycs  =  3.6,  the  rate  ratios  shown 
in  the  last  column  are  obtained  as  outlined  above.  Clearly,  as  the  AJ  =  3 
channel  becomes  open  (at  J=31),  f  is  small  and  the  corrections  to  the 
RA  model  are  minor.  As  f  increases  (i.e.,  as  the  AJ=3  channel  opens 
more  fully),  the  corrections  grow.  The  qualitatively  important  feature 
suggested  by  these  results  is  that  the  autoionization  rates  should 
increase  with  J  more  rapidly  than  predicted  in  the  original  RA  model 
whenever  new  channel  openings  occur.  In  Figure  2  below,  the 
linewidths  estimated  by  Clary  for  these  transitions  and  those 
obtained  after  applying  the  above  frustrated-channel  corrections  are 
shown.  The  quantitive  predictions  achieved  through  these  corrections 
probably  become  less  reliable  as  the  corrections  become  larger 
because  the  perturbative  method  used  to  estimate  such  corrections 
becomes  susppect. 

VII.  Summary 

Modifications  have  been  made  to  Clary’s  original  rotationally 
adiabatic  model  that  treat  final-state  product  yields  and  rate 
enhancements  near  channel  openings.  The  latter  are  found,  in  an 
example  calculation  chosen  to  represent  rotational  autoionization  of 
dipole-bound  CH2CN-  (with  EA=60  cm-l  and  PDebyes=3.6)  to  produce 
substantial  rate  enhancements  in  the  neighborhood  of  new  channel 
openings.  The  essential  ingredients  in  this  model  are  the  electron 
binding  energy,  the  dipole  moment,  and  the  rotational  constants  of 
the  molecule.  The  corrections  are  tabulated  here  in  a  manner  that 
facilitates  application  to  a  wide  variety  of  dipole-bound  anions.  For 
any  anion  state  J,  if  the  AJ  value  to  the  first  energetically  accessible 
neutral-molecule  state  and  the  anion-neutral  state  energy  difference 
are  knov/n,  then  T(n,f)  and  W(n,J)  can  easily  be  calculated.  The  tables 
provided  here  then  permit  easy  evaluation  of  the  corrections  to  the 
rotationally  adiabatic  rates. 
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2.  An  analogous  treatment  for  symmetric  top  molecules  coupled  to  a 
dipole-bound  electron  is  contained  in  Clary's  original  works for 
pedagogical  purposes,  emphasis  is  placed  here  on  the  diatomic- 
molecule  case  where  the  notation  is  less  complicated.  The 
enhancements  to  the  original  Clary  model  introduced  here  for  the 
diatomic-molecule  case  apply  equally  well  to  the  symmetric  top  and 
other  situations. 

3.  Smaller  contributions  from  j=J-2,l=2  ,  j=J-3,l=3,  and  so  on  are,  of 
course,  also  present  i.i  this  hybrid  at  these  R-values, 

4.  For  small  n,  this  places  the  crossing  between  the  the  neutral's  (J-n) 
and  (J-n-i-1)  states;  however,  for  n  >  3  the  crossing  can  occur  above 
even  the  (J-n+1)  state. 

5.  The  experimentally  observed  small  differences  in  the  rotational 
constants  of  the  neutral  and  the  dipole-bound  anion,  caused  by 
geometrical  differences  induced  by  the  "extra"  electron’s  bonding  or 
antibonding  effects,  give  rise  to  progressions  of  rotational  energy 
levels  which  are  not  adequately  described  within  the  above  picture. 
For  example,  if  the  neutral  molecule  has  a  by  value  (b^)  which  is 
slightly  larger  than  that  of  the  dipole-bound  anion  (b’)  and  this  state 
of  the  anion  is  bound  by  EA,  then  the  pattern  of  neutral  and  anion 
states  with  energies 

EOj  =  bO  J(J-Hl)  and 
E-j-  =  b-  J'(J'+1)  -  EA 

gives  rise  to  a  situation  in  which  the  lowest  states  of  the  anion  that 
can  detach  can  involve  relatively  large  AJ  .  As  J'  is  increased,  the 
values  of  AJ  then  decrease  (  because  the  neutral's  rotational  energy 
level  ladder  moves  faster  than  does  the  anion's).  However,  as  one 
moves  to  higher  and  higher  J’,  there  is  a  minimum  value  (  whose 
value  depends  on  EA,  b',  and  b^)  which  AJ  achieves;  after  this 
minimum  value,  AJ  increases  with  increasing  J'.  Thus,  unlike  the 
bO=b  *  case,  one  need  not  realize  AJ  =  1  at  all.;  if  EA  is  large,  the 
minimum  AJ  can  be  quite  large,  in  which  case  detachment  will  be 
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slow  for  all  J'  states.  In  contrast,  if  b*  is  slightly  larger  than  b®  ,  the 
lowest  states  of  the  dipole-bound  anion  that  can  detach  undergo 
larger  AJ  transitions.  As  J’  increases,  smaller  and  smaller  AJ  values 
are  realized.  Eventually  AJ  =  0  transitions  open,  and  at  even  higher  J' 
values,  AJ  becomes  negative  (  i.e.,  transitions  in  which  the  neutral  is 
in  a  higher  rotational  level  than  the  anion  occur)  and  becomes  even 
more  negative  as  J'  further  increases.  Unlike  the  bO  =  b-  case,  AJ  =0 
and  even  negative  AJ  values  are  realized  here. 

6.  L.  D.  Landau,  Phys.  Z.  Sow.,  1,  46  (1932);  C.  Zener,  Proc.  Roy.  Soc. 

A137.  696  (1932);  Stueckelberg,  E.  C.  G.,  Helv.  Phys.  Acta.  5.,  369 
(1932). 


Figure  1.  Rotationally  diabatic  and  adiabatic  potentials  for  neutral- 
molecule  rotor  and  electronic  angular  momenta  j,l  =  34,0;  33,1;  32,2; 
31,3;  30,4;  29,5;  28,6;  and  27,7  for  CH2CN-  with  K=0.  The  anion  level 
J=34  is  also  shown.  See  text  for  details. 

Figure  2.  Linewidths  in  Mhz  for  CH2CN-  for  J=31  to  39  within  Clary's 
adiabatic  model  and  the  present  corrected  model. 


Table  1 

Anion  and  Neutral  Rotational  Energies  (cm'^)  for  CH2CN‘  ,  K=0; 

b=  0.341  cm-1,  EA=60  cm-^ 


Anion 

Neutral 

J 

AJa 

-60 

0 

0 

- 

-58 

2 

2 

- 

-53 

7 

4 

- 

-46 

1  4 

6 

- 

-35 

24 

8 

- 

-22 

38 

10 

- 

-7 

53 

1  2 

- 

12 

72 

14 

5 

33 

93 

16 

5 

57 

117 

18 

5 

83 

143 

20 

5 

113 

173 

22 

5 

145 

205 

24 

4 

179 

239 

26 

4 

217 

277 

28 

4 

257 

317 

30 

4 

301 

360 

32 

3 

346 

406 

34 

3 

394 

454 

36 

3 

445 

505 

38 

3 

a.  AJ  is  the  change  in  J  accompanying  a  transition 

from  this  J-level 

the  anion  to 

the  first  energetically 

accessible  state 

of  the  neutral. 

Table  2 

T(n,0  Factors  for  AJ=n  Values  from  1  to  4  and  f-Values  from  0.1  to  1 


T(n,f)=(l-f)2  Xn(f) 


f= 

n=l 

n=2 

n=3 

n=4 

0.1 

2.4x10-3 

3.5x10-5 

5.5x10-8 

9x10-9 

0.2 

0.027 

3.8x10-3 

3x10-4 

3x10-5 

0.3 

0.070 

0.017 

0.004 

0.001 

0.4 

0.100 

0.040 

0.016 

0.006 

0.5 

0.110 

0.060 

0.034 

0.019 

0.6 

0.096 

0.066 

0.046 

0.032 

0.7 

0.066 

0.053 

0.042 

0.034 

0.8 

0.035 

0.031 

0.028 

0.026 

0.9 

0.010 

0.009 

0.009 

0.008 

Table  3 

W(n,J)/|i2  Factors  for  AJ=n  from  1  to  4 

n=l  32  (J-0.5)/[n2Debyes(J-l)l 

n=2  271  (J-1.5)/[ti2Dcbycs(J*2)] 

n=3  1047  (J-2.5)/[^2Debycs(J-3)] 

n=4  2845  (J-3.5)/[^i2Debyes(J-4)] 
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Table  4 

Anion&Neutral  Energies  (cm'^),  f-Values,  and  Rate  Ratios  for^  CH2CN' 


Anion 

Anion 

Neutral 

f- 

Rate 

J  value 

Energy 

Energy 

Value 

Ratioh 

3  1 

278.27 

276.89 

0.05 

1.00 

32 

300.10 

296.67 

0.13 

1.00 

33 

322.60 

317.13 

0.18 

1.02 

34 

345.79 

338.27 

0.24 

1.16 

35 

369.66 

360.10 

0.29 

1.31 

36 

394.21 

382.60 

0.34 

1.78 

37 

419.45 

405.79 

0.39 

2.25 

38 

445.36 

429.66 

0.44 

2.96 

39 

471.96 

454.21 

0.48 

3.57 

40 

499.24 

479.45 

0.52 

4.04 

a.  A  pseudo-diatomic  b 

value  of  0.341  cm-l 

and  an  EA  of 

60  cm'l 

were  used  to 

compute 

these  energy 

levels.  They 

correspond  to  K'=0  =>  K=0  transitions  of  the  nearly 
symmetric-top  CH2CN*  species.  For  all  anion-to-neutral 
transitions,  AJ=3,  so  the  neutral  J  is  equal  to  the  tabulated 
anion-J  minus  3. 

b.  Defined  as  l+T(n,f)W(n,J)/p2Dg{jyg5^  this  is  the  "corrected" - 
to-rotationally-adiabatic  rate  ratio. 
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